Exact interferometers for the concurrence and residual 3-tangle 
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In this paper we describe a set of circuits that can measure the concurrence of a two qubit density 
matrix without requiring the deliberate addition of noise. We then extend these methods to obtain 
a circuit to measure one type of three qubit entanglement for pure states, namely the 3-tangle. 

PACS numbers: 03.67.Mn, 03.67.-a, 03.65.-w 



Introduction Techniques for measuring entanglement 
without first reconstructing the state have recently at- 
tracted a lot of attention. Some of these approaches 
have been based on the Structural Physical Approxima- 
tion (SPA) [if, followed by measuring the spectrum of the 
resulting density matrix 0,0,01 (but see also d,|(J). One 
method for measuring that spectrum was inspired by 
in which it was shown that while the physical evolution 
of the density matrix must be p h> U pU\ the interfer- 
ence pattern is proportional to Re(ve~' lip ) =Re(Tr(C/p)), 
where v is the visibility and tp is the phase shift. Thus the 
interferometer circuit in Figure n can modify the expec- 
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FIG. 1: General interferometer circuit 

tation value of the measured qubit, which corresponds to 
the path in a Mach-Zehnder interferometer. Note that U 
can be any unitary operation. 

In |j| we showed how to measure the spectrum of p T2 
without using the SPA to shift the spectrum to be non- 
negative. The "Wootters spin-flip" used to define the 
concurrence is also an unphysical operation. The tech- 
niques in cannot be used directly to measure this in- 
variant, but a more general type of circuit can. 

We can measure the 3-tangle in a similar way, but 
some of the algebra used in the design of these circuits is 
only valid for pure states; the function conflates subsys- 
tem mixing with entanglement, and so these circuits will 
have the same restrictions on their validity. 

A Generalized Interferometer We will begin by gen- 
eralizing the circuit construction in 0. This studied 
the behaviour of a Mach-Zehnder interferometer. By as- 
sumption, our particle has an internal degree of freedom, 
or spin. The circuit couples the internal and external 



degrees of freedom, by the operation 



e lx 




(1) 



When the control qubit is "on" , the internal degrees of 
freedom p undergo the evolution p t— * U^'PqU^ where 
the label (i) signifies that only the internal degrees of 
freedom are affected. Note that can be any unitary 
matrix that only acts on the internal degrees of freedom. 
A little algebra then gives us that the output intensity, 

I oc 1 + |Tr([/»p )| cos[ X - argTr([/«p )]. (2) 

We will now modify this circuit to generalise Q to a 
completely positive (CP) map. (While this question has 
been considered before in [lOj , this had a different moti- 
vation and was not concerned with measuring entangle- 
ment.) We will limit ourselves to convex combinations 
of unitary operators, and introduce a family of unitary 
transformations labelled by the index fc, 



U k = 




1 



® u, 



e lx 




1 



(3) 



Note that these operators differ only in the unitary op- 
erator U k acting on the internal degrees of freedom. If 
we weight these with probabilities pk , then the controlled 
evolution of the internal degree of freedom is now 



k 



(4) 



When we redo the calculation in with our general- 
ized internal dynamics, we obtain the measured intensity 
along the |0)(0| arm of the interferometer, which will be 
proportional to 



argTr 5> fe [/£ 



Pa 



(5) 

We will now use this to measure the concurrence. 

Measuring the Concurrence The concurrence is a 
bipartite entanglement measure defined on mixed states 
of two qubits to be 



Cab = max{Ai - A 2 - A 3 - A 4 , 0} 



(6) 
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where the AjS are, in decreasing order of magnitude, the 
positive square roots of the eigenvalues of pab Pab, where 
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(7) 



and where T denotes matrix transposition in the compu- 
tational basis. The method in || will not work, because 
we will obtain acausal circuits in which some of the in- 
put rails double back on each other. If we had wanted 
the spectrum of pab on its own, we could measure that 
in the same way as in t 8|. However, when we try to 
pre-multiply by pab, we find that we have to sum two 
pairs of input indices with two other input indices and 
likewise with two pairs of output indices. The SPA ^1 
would also have serious problems because of that matrix 
multiplication, as it will instead measure the spectrum of 



PabHpab) = PPab + vpAB<J y ® o y p\ B a y ® a y , (8) 

where p >> v arc the parameters defining the SPA for 
this map. In effect, the SPA circuits will measure the 
eigenspectrum of pab , with a small correction term pro- 
portional to pabPab- Without prior knowledge of pab 
there is no way to correct this error. 

We will now describe two alternative methods for mea- 
suring the concurrence. Recall that 
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= ie, where e = 
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(9) 



Rewriting equation Q in index notation and using the 
Einstein summation convention, we obtain 

(% ={e®e)^{p T f r l{e®e)f q , (10) 

where we have omitted the AB subscripts to reduce clut- 
ter. Converting the e matrices to vectors, we can write 
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Now use the identity e ab e cd = 5%5$ - 5%5% to obtain |l3| 

/-!;; i' 1 ;: i>7/; ^;/C • m 

In the more familiar matrix notation, this is 

Pab — Pab ~ ^a ® Pb ~ Pa <£> Is + 1ab< (13) 

We can now safely perform the matrix multiplication: 

PabPab = Pab ~ Pab(^a® Pb) ~ Pab(pa®^b) + Pab- 

(14) 

We will need to find a way to measure the various mo- 
ments of pabPab to determine the eigenspectrum and 
thus obtain the A^s for ©. There are two ways in which 
we can convert equation l|14f> into an interferometer. The 
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FIG. 2: Summation diagram for pabPab- 



first method uses the modified interferometry circuit in- 
troduced in the previous section directly, as it enables us 
to measure the effects of sums of unitary matrices, de- 
spite the fact that no physical evolution can be written 
as the sum of unitary matrices in this way. It should be 
noted that we will be needing the moments of pp, rather 
than p itself. So, the fact that some of the terms in p 
do not have trace 1 will not require us to rescale those 
terms, as p\ — p, so the second set of rails can just be 
ignored, without needing to find a way to change the 
relative weights of the terms. 

The permutations required for each term can be de- 
composed into two operations which are variously "on" 
or "off" for each term in l|14|) . The set of probabilities 
also factorizes conveniently, so we will need only two aux- 
illiary control ancillas. Figure [5] shows the terms of i|14|) 
as a summation diagram, which is a pictorial version of 
index notation. The nodes represent indices and the links 
represent summations. Unlinked nodes are free indices. 
The nodes on each copy of p are ordered as for standard 
tensor index notation. The summation pattern for the 
last term is defined to have both ancillas in the "off" po- 
sition. The corresponding circuit diagram for this single 
term could be obtained by rotating the Figure clockwise 
by —90°, omitting the p symbols and inserting the result- 
ing wiring pattern into the box marked "U" in Figure ^ 
In this case, all the rails go straight through and this sub- 
circuit would measure the trace norm if implemented on 
its own. Contrast this with the first term which will have 
both ancillas "on" and is proportional to Tr(p AB ). Note 
that the only difference between these two terms is a pair 
of swaps. 

We must also take care to ensure that we reproduce 
the minus signs in equation Q14[l. so the qubit we are 
going to measure will be able to "see" the minus signs 
from the ancillas. This can be done either by initializing 
the control ancillas with the appropiate relative phase, 
together with an initial controlled-dz gate, or by includ- 
ing the desired phase in the gate itself We will use 
the first method and we will therefore need to initialize 
both auxilliary ancillas in the state — ) = ^(|0) — |1)). 
The fact that the ancillas must be normalized means that 
each auxiliary ancilla will decrease the visibility of the in- 
terference fringes by a factor of 2. If we now insert these 
components into a modified version of Figure ^ so they 
are controlled by the two auxilliary ancillas, we will ob- 
tain the full circuit for Tr (pab Pab) m Figure |3| 

The same principles can be used to construct the cir- 
cuit for Tr ((pab Pab) 2 )- This has 16 terms and will re- 
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FIG. 3: The circuit for jjTrpAB/311?. 



quire four auxilliary control ancillas, initialized in the 
same way. The structure of this circuit can be parsed 
as follows. The box in dotted lines with no auxilliary 
controls performs the permutation for the matrix multi- 
plication of the two copies of pabPab- The term with all 
the auxilliary control ancillas off will look like Figure 0] 



FIG. 4: All auxilliary control ancillas off. 

Figure [S] is the diagram for the term with all the ancilla- 
controlled operations "on". Again, the difference be- 
tween the two is just a set of simple swaps, 



FIG. 5: All auxilliary control ancillas on. 

which can be controlled by auxilliary ancillas as before. 
Thus we obtain the circuit for measuring ^Tr((pp) 2 ), 
in Figure The circuits for ^Tr ((pabPab) 3 ) and 
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FIG. 6: The circuit for ^Tr^pp) 2 ). 

i5^Tr((pAB/9As) 4 ) can be constructed in a similar way. 
We can now calculate the spectrum {r]i} of pabPab us- 
ing the recipe given in 0, 0] . The positive square roots 



of the 7/iS give the concurrence spectrum and equation 
© can then be used to compute the concurrence. 

This first method shows that the concurrence (and 
hence the entanglement of formation) can be measured 
without any prior knowledge of the state. However, this 
method is not very sensitive, because the visibility drops 
by a factor of i for the mth order concurrence moment. 
Fortunately, there is a second method for measuring the 
moments of the concurrence spectrum. This uses the 
principle of implicit measurement [lfij . whereby we can 
assume the auxilliary ancillas have been measured at the 
end of the circuit without affecting the statistics of the 
principal measurement. These circuits are examples of 
"non-erasing quantum erasers" |l7j . albiet very simple 
ones. In effect we will be measuring each term in the 
various summations separately. If we choose to measure 
the ancillas, we will be partitioning the output of the 
primary ancilla into subensembles with greatly enhanced 
visibility. We can also simplify the circuits by initializ- 
ing the auxilliary ancillas using only a simple Hadamard 
gate, and omitting the controlled-cz gates altogether. 

Alternatively, we can implement the circuits for each 
term separately, and insert the necessary weight factors 
by hand. This has the additional advantage that we 
can exploit the cyclicity of the trace and the fact that 
^Pab = 1 and its first four moments completely define 
its spectrum, to reduce the number of circuits required 
from 340 to no more than 111. This does not compare well 
to full state tomography (which needs only 16 parame- 
ters) but this upper bound is unlikely to be tight, as this 
system can be characterised by 18 non-local parameters 
|l8j . of which only 9 are fully algebraically independent; 
the other 9 are discrete valued. 

Measuring the 3-tangle The residual 3-tangle as de- 
fined in [9j is the modulus of a complex number. It is only 
defined in closed form for pure states of three qubits. If 
we define the 2-tangle tab to be the square of the con- 
currence 0, then the residual 3-tangle is defined as en- 
tanglement between the qubits that cannot be accounted 
for in terms of bipartite entanglement: 



TABC = Ta(BC) — T AB ~ TAG- 



(15) 



The measure is symmetric under relabelling of the three 
parties 0|. We cannot use (fT5)l directly, because we can- 
not define the equivalent of a y on the merged party (BC) 
without knowing the support of PbCi which would require 
detailed knowledge of the state. 

The fact that the two-party density matrices are all 
rank 2 means that A3 = A4 = 0. So we can write 



tag = (Ai - A 2 ) 2 = A 2 - A 2 , - 2AiA 2 . 



(16) 



Wootters et al. then obtain that ta(bc) = 4det pa- This, 
together with the definition in l|15fl gives us that 



tab + tac = 4 det p A ~ t A bc- 



(17) 
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Using the relabelling symmetry and the fact that detpj = 
5(1 — Tr(p?)) for 2 x 2 matrices, we can write 

t abc = 2(1 - Tr(p 2 A ) - Tr(p|) + Tr(p* ) - t ab ). (18) 

We can calculate this from the rescaled expectation val- 
ues measured by the circuits for tab we obtained ear- 
lier, together with the statistics from the circuits for the 
Tr(p|)s. These one-party circuits need two copies of their 
single particle ps each and the controlled- U operation for 
these is the SWAP gate [2J . We can combine these three 
one-party circuits into a single three party circuit with 
three control qubits, thus saving four copies of P123 per 
measurement. We can therefore measure the 3-tangle us- 
ing this family of five circuits, provided the state is pure. 

We can improve both the cost and sensitivity for the 
3-tangle as follows. Recall that [j| 

Tt(pabPab) + Tr{pAcPAC') = 4detpA, (19) 

where the superscripts refer to which two-party density 
matrix was used to obtain the concurrence spectrum. We 
know that \[ AB) \^ B) = \ {BC) X { 2 BC) = \[ AC) \i AC) from 
the relabelling symmetry, so we can omit the superscripts 
and write tabc = 4A1A2. The square root in the defini- 
tion of the AiS prevents us from measuring this but we 
can construct a circuit for |tabc| 2 = lOA^A^. Now we 
just need to notice that 2Af A| = (A 2 + A 2 .) 2 — Af - A|, 
and then we can write 

\tabc\ 2 = 8 (^{pabpTb? - Tt(( P abpIb) 2 )) ■ (20) 

A circuit with visibility oc | tabc 1 2 can be obtained from 
Figure HO We need an equal superposition of the sub- 
circuits with the matrix multiplication (the faint dotted 
box) switched on and off so we will insert a new ancilla 
in the state |— ) and another controlled-o^ as we want 
the term with the matrix multiplication on to acquire a 
minus sign. The intensity along the |0)(0| arm for this 
naive circuit will be |jj 

I 1tabc1 2<xI + \tabc\ 2 /32. (21) 

This circuit needs four copies of the state per run, com- 
pared with the first method which needs 22 copies. The 
sensitivity of this naive circuit can also be improved by 
measuring the terms separately. This can in fact be done 
with just 14 circuits, compared with the 16 parameters 
that must be determined to learn the state of a three 
qubit state that is known to be pure 20] . 

Discussion We have presented circuits for both the 
concurrence and the residual 3-tangle. The circuits for 
the concurrence should work for any two qubit state. The 
visibility of these circuits is not much better than that for 
the SPA circuits, however, the former do not introduce 
state-dependent errors. The visibility can be improved 
by decomposing the original moment circuits into unitary 
circuits, but there will then be many more of them. 



The two methods for finding the residual 3-tangle are, 
strictly speaking, only valid for pure states. In any lab- 
oratory setting the states will have a certain amount of 
mixing. The effect of this on the residual 3-tangle is in- 
completely understood so it is unclear how useful these 
circuits would be in practice, but at least these circuits do 
not aggravate the mixing beyond that caused by unavoid- 
able experimental noise. They also require fewer circuits 
than the corresponding state tomography problem. 

We will generalize the concurrence circuits to measure 
other multipartite entanglement monotones in further 
work j2l|]. 
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